Charged Particle Ratio Fluctuation as a Signal for QGP 



S. Jeon* and V. Koch^ 
Nuclear Science Division 
Lawrence Berkeley National Laboratory 
Berkeley, CA 94720, USA 

In this letter we argue that the event-by-event fluctuations of the ratio of the positively charged 
and the negatively charged pions provides a signal of quark-gluon plasma. The fact that quarks 
carry fractional charges is ultimately responsible for this distinct signal. 
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It is of great importance that we have a clear signal of 
the long-sought quark-gluon plasma (QGP) not only for 
the experiments at RHIC but also for theoretical reasons. 
At stake is our fundamental understanding of strong in- 
teractions as well as understanding of the state of matter 
in the very early universe fjl]. Proposed signals of this 
new state of matter abound in literature 2] one of the 
most studied being the J/tp suppression Sfl. 

In this paper, we propose the event-by-event h + /h~ 
fluctuations as a distinct signal of QGP formation. We 
would also like to stress that this observable is something 
that can be and already has been calculated on a lattice. 

The idea is very simple and is reminiscent of the origi- 
nal detection of color in e + e~ experiment where one mea- 
sures 
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Hadrons 
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(i) 



Here Q q is the charge of each flavor and iV c is the number 
of colors. Note that if the fundamental degrees of free- 
dom were hadrons, R e + e - would be very different from 
this simple counting. We would like to establish that 
the event-by-event h + /h~ fluctuations can similarly de- 
termine whether the underlying degrees of freedom are 
quarks and gluons or hadrons. 

The point is that in the QGP phase, the unit of charge 
is 1/3 while in the hadronic phase, the unit of charge 
is 1. The net charge, of course does not depend on such 
subtleties. However, the fluctuation in the net charge de- 
pends on the squares of the charges and hence strongly 
depend on which phase it originates from. Measuring the 
charge fluctuation itself, however, is plagued by system- 
atic uncertainties such as volume fluctuations due to the 
impact parameter variation. In a previous letter |5|, we 
showed that the multiplicity ratio fluctuation is only sen- 
sitive to the density fluctuations and not to the volume 
fluctuations. The task for us is then to find a suitable 
ratio whose fluctuation is easy to measure and simply 
related to the net charge fluctuation. 

The obvious candidate is the ratio F = Q/N c b where 



Q = N + - N- 
is the net charge and 



N, 



ch 



N+ + N_ 



(2) 



(3) 



is the charge multiplicity. Here N± denote the positive 
and negative multiplicities. Instead of using F, how- 
ever, in this paper we propose to use the charge ratio 
R = N+/N- . The advantages of using R over F are that 
although trivially related, R is more fundamental to ex- 
periments and the signal is about 4 times amplified in R 
as we show below. 

To relate R with F, we first rewrite the charge ratio as 



R = 



1+F 

1 - F 



(4) 



When (Neb) 3> (Q) we can safely say |F| <C 1. Expand- 
ing in terms of F yields 



R « 1 + 2F + 2F 2 



(5) 



Defining Sx = x — (x) for any fluctuating quantity x, it 
is easy to show 



(SR 2 ) = (R 2 )-(R) 2 ^4(SF 2 ) 



(6) 



where (• • •} denotes the average over all events. 

Let us now consider (SF 2 ^ more closely. In a previous 
letter || (see also |§ and the upcoming paper 0), we 
showed that a ratio fluctuation can be expressed as 



(SF 2 ) 



5N, 



ch 



(Q) 

(N ch f \\(Q) (N) ch/ 



(7) 



We then showed that when the average ratio is very much 
different than 1, the fluctuation is driven mainly by the 
fluctuation in the smaller quantity (for instance K/~k fluc- 
tuation is driven by K fluctuation). At RHIC we expect 
(Q) I (N c b) ~ 5%. Hence the fluctuation in F is totally 
dominated by the fluctuation in Q so that 
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where N q+q from here on denotes the number of quarks 
and anti-quarks. The constant 



If we can detect all charged particles from a heavy ion 
collision, the net charge Q is a fixed quantity and hence 
will not fluctuate. This implies that (SF 2 ) is very small 
with a 47r coverage. However, no detector can catch all 
charged particles. Our study dj7j] shows that for a real- 
istic detector acceptance, using the grand canonical en- 
semble is acceptable and that is what we assume here. 
The corrections to this approximation have been worked 
out and will be reported in the upcoming paper R] . (Also 
see Finite Acceptance Correction at the end of this 
letter.) Our main observable is then, 

(N ch ) (SR 2 ) = 4<A ch > (5F*) = 4^ (9) 

to the leading order in the fluctuations and 1/ (Nch). 

So far, we have only considered statistics of the ratio 
fluctuations. Physics lies in how the charge fluctuation is 
expressed in terms of the fluctuations in the fundamental 
degrees of freedom. For simplicity, let us consider a pion 
gas and a QGP consisting of u and d quarks and gluons. 
Our main conclusion does not depend on this simplify- 
ing assumption. We will briefly consider the size of the 
corrections towards the end of the paper. 

In a pion gas, the fundamental degrees of freedom are 
of course pions. Hence, A c h = N TT + + N n - and 

5Q = <5A 7r+ - SN^- (10) 

Using thermal distributions and disregarding correla- 
tions, we get 

(5Q 2 ) = (SNl) + (SN 2 ) = w„ (N ch ) (11) 

where 

w v ee (5N$) I (N„) (12) 

is slightly bigger than 1 J|||. Hence for a pion gas, 

D had EE (7V ch )(^ 2 )| hadron «4. (13) 

For a quark-gluon plasma, 

SQ = Q u S (N u - N u ) + Q d S (N d - N s ) (14) 

where Q q is the charges of the quarks and N q is the num- 
ber of quarks. The fluctuations (SN^ are measured on 

lattice H and we will shortly get back to the results. For 
now, let us consider a thermalized gas of non-interacting 
quarks and gluons to get a physical baseline. Thermal 
distributions and no correlations yield 

(SQ 2 ) = Qlw u (N u+a ) + Q 2 d w d (N d+d ) (15) 



w q ee (SN 2 ) I (N g ) (16) 

is slightly smaller than 1 due to the fcrmionic nature of 
the quarks. 

Relating the final charged particle multiplicity A c h to 
the number of primordial quarks and gluons is not as 
simple. To make an estimate, we assume that the en- 
tropy is conserved || and that all the particles involved 
are massless, in thermal equilibrium and non-interacting. 
For such particles, the following relation between the en- 
tropy density and the particle number density holds: 

a B = 3.6n B (17) 

and 

cr F =4.2n F (18) 

where the subscript B, F signifies the particle types. The 
total entropy of a quark-gluon gas in a given volume V qg 
is 

S = V qg a qg = 3.6 (N g ) + 4.2 ((N u+B ) + (N d+d )) (19) 

where N g is the number of gluons inside the volume and 
N q + q is the number of quarks and anti-quarks inside the 
same volume. As the volume expands and cools, eventu- 
ally the quarks and gluons are converted to pions. Since 
entropy is conserved, the number of pions coming from 
these quarks and gluons must be given by 

W = ^ = (N g ) + || ((N u +u) + (N d+d )) (20) 

The charged multiplicity is 2/3 of N^ due to isospin sym- 
metry: 

(N c h) = \ {(N g ) + 1.2 (N u+a ) + 1.2 (N d+d )) (21) 

Then for massless non-interacting quarks and gluons, 

-Dqgp = (A ch ) <<5i? 2 )| QGP « 0.75 (22) 

from Eq.(||) and using Eqs.(||) and (|l]). The value of 
-Dqgp is more than a factor of 5 smaller than the value of 
Z?had in Eq. (|l3|) ! This is an unmistakable signal of QGP 
formation from such a simple measurement. 

We now would like to stress that (5Q 2 ) / (N c h) is al- 
ready calculated on lattice and hence one does not have 
to rely on the above thermal model calculation. In Ref. 

, Gottlieb et al report their calculation of the quark 
number susceptibility and the entropy density with 2 fla- 
vors of dynamic quarks. These two quantities are directly 
related to the net charge fluctuation and the charged mul- 
tiplicity in the following way. 

From the definition of the charge susceptibility Xqi it 
is clear that 
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(SQ 2 )=V qg T Xq 



(23) 



where V qg is the volume of the quark-gluon plasma at the 
hadronization and T is the temperature. Gottlieb et al 
calculated the quark number density susceptibilities 



T\s = V qg l(5n u + Sn d y 



and 



Txns = V qg l(5n u - 5n d ) 



(24) 
(25) 



and found that at high temperature both are very close 
to the non-interacting thermal gas limit 



Xs ~ Xns 



IT 1 



(26) 



From this result, one can first of all infer that u and d 
quark densities are uncorrelated 



and 



(SN u 6N d ) « 
(8N 2 ) » (8NI) 



(27) 
(28) 



These results imply that the charge fluctuation at high 
temperature follows that of the thermal fermion gas 



or 



(29) 
(30) 



For the charged multiplicity, we assume that the rela- 
tions Eqs. (||^ - ^l|) still hold. Equating the entropy of 
the final pions with that of the primordial quarks and 
gluons, one obtains 



(iVch) 



1 

5A 



Vqg (cr qg ) 



(31) 



Ref. H reports that the gluon entropy density <r g is al- 
most the same as the non-interacting thermal bosons, but 
the quark entropy density a u+d is about one half of that 
of the non-interacting thermal fermions. Hence, the total 
entropy from the lattice calculation is 

{(Tqg) ={(Tg)+ ((Tu+d) 

= 16 x 3.6 x (f g ) + 24 x 4.2 x a (f q ) 
«12T 3 (32) 

where a w 1/2 and (f q . g ) is the average density per de- 
gree of freedom. 

Using Eqs.(p0| - B2), the lattice calculation gives 



Aat = (N ch ) (5R 2 ) 



lattice 



l 



(33) 



which is still 4 times smaller than the pion gas result. 
This is the main result of this paper. The difference be- 
tween the pion gas result and the QGP result is distinct 
enough one should easily see it in the first few days of 



data collecting at RHIC. We are of course aware that lat- 
tice result is not yet exact. We hope that this discussion 
will actually stimulate more sophisticated lattice calcu- 
lations of quark number susceptibilities. We also note 
that this is an opportunity to test what is calculated on 
lattice in a direct observation. 

The picture obtained above holds if the following 
two conditions are met: (i) The detected phase space 
is a small sub-system of the whole, (ii) The original 
quarks and gluons stay in the system during or after the 
hadronization. Both conditions can be met if the rapidity 
intervals are such that 



2/total ^* 2/accept > 1 



(34) 



Here, ytotai is the rapidity range allowed by the energetics 
of the collisions and y accC pt is the acceptance interval of 
a given detector. The first of these conditions is needed 
to ensure that the rest of the system acts as a reservoir 
and the second condition ensures that the charge diffu- 
sion in the rapidity space during and after hadronization 
is negligible. In real life, of course, Eq. ([mJ) is satisfied 
in varying degrees. For instance, the STAR at RHIC has 
ytotai ~ 10, 2/acccpt ~ 3 and hence corrections should be 
taken into account. We would like now to discuss Caveats 
and corrections due to these and other effects. 
Hadronization : First, charge conservation fixes the 
net charge once it is set in the QGP phase. Subsequent 
hadronization cannot change the net charge. Second, 
the entropy can only increase during the hadronization. 
Hence, our estimate of (SQ 2 ) / (Neb) from QGP should 
not only survive the hadronization but also could even 
become smaller, thus strengthening the signal. 
Phase Space Cuts : A possible issue of having a rapid- 
ity cut is that this implies a box in the momentum space 
while the argument presented so far dealt with a box in 
the coordinate space. Within the Bjorken scenario, this 
is of course not a problem. However this is not a problem 
even without the Bjorken picture. If anything, it actu- 
ally helps. Once it is established that the ratio h + /h~ 
is independent of the overall volume, the only crucial in- 
gredient in our argument is that the number fluctuations 
are all Poissonian, namely, 



(SN 2 )/(N) 



1 



The quantum statistics makes this ratio deviate from 1. 
For bosons (pions), restricting to small momenta or ra- 
pidity makes this ratio bigger. For fermions (quarks), it 
makes this ratio smaller. From Eqs.([TT|), (fl^), (15) and 
( |l6| ) , it is clear that having a momentum cut-off can then 
only enhance the contrast. 

Resonance Contributions : As explained in a previ- 
ous paper H , neutral resonances introduce positive cor- 
relations between N + and iV_ and hence lower the value 
of -Dhad from 4. In a thermal scenario studied in the same 
paper, we found that the resonances reduce the fluctua- 
tion by about 30 %. Hence for a realistic hadron gas, 



3 



D 



had 



(N ch ) (5R 2 ) 



(35) 



This is still a factor of 3 bigger than the lattice result. 
Mixtures : If the system is a mixture of a QGP and a 
hadron gas, the signal should depend on the fractions. To 
a first approximation, it should be a linear combination 
of Eqs. @ and © 

<iV ch ) (SR 2 ) = D had (l - /) + Aat/ 

where / is the QGP fraction. Even if / ~ 0.5, the signal 
can be still visible. 

Rapidity Correlations : If gluon interactions dom- 
inate the creation of hadrons in high energy hadron- 
hadron collisions, the unlike charged particles are 
strongly correlated in rapidity. Such strong correlation 
can further lower the value of -Dhad- If the reduction is 
big enough to mimic QGP signal, we should certainly see 
the false signal in even in the nucleon-nucleon collisions. 
Using the two particle correlation data in pp collisions 
compiled by Whitmore JlQ] , we performed Monte Carlo 
simulations for such a scenario. The correction to ([35]) is 
found to be small (~ 10 %). This is due to the experimen- 
tal fact that even though the correlation is pronounced 
in the 'connected' part of the correlation i?(y + ,j/_), it 
is not so pronounced in the full two particle correlation 
function 



P2(y+,y-) = pi(y+)pi(y-) [l + R(y+,y-)} 



(36) 



For details, see pj. 

Finite Acceptance Correction : The finite size of the 
acceptance window introduces a factor of (1 — p) cor- 
rections where p is the fraction of the total multiplicity 
inside the acceptance window. This is easy to under- 
stand. If the detector sees a 100 % of all charged parti- 
cles {p = 1), the fluctuations should shrink down to zero 
due to the global charge conservation. Fortunately, this 
is a common factor that applies to both Eqs. ( |35f ) and 
( |33| ) so that the ratio stays the same. More details will 
be reported in 0. 

Effects of Rescattering : The partons as well as 
hadrons are subject to rescattering during the course of 
a heavy ion collision. This in principle may affect the 
above charge fluctuations by diffusing the charge in ra- 
pidity space. However in the limit y total > J/accopt > 1 
these effects should be very small since they scale as the 
surface to volume ratio in rapidity space. To estimate 
the effect in the hadronic phase, we performed a sim- 
ple Monte Carlo calculation where a Gaussian noise with 
a = 0.5 are added to each particle's rapidity originating 
from a highly correlated source. Doing so increases the 
value of (N ch ) (SR 2 ) up to 40 % assuming the STAR ac- 
ceptance at RHIC. In a slightly different context, Gavin 
jnj also found that when a strong flow is established, the 
rescattering effect on the hadron rapidities is relatively 
small. 



Strangeness : Adding Kaons to a pion gas will not 
change the value of £>h a d because their contribution 
adds exactly the same amount to both numerator (c.f. 
Eq. (pi]) ) and the denominator. For a quark gluon plasma, 
the lattice calculation |)| suggests that at high tempera- 
ture, the strangeness entropy is about 40 % of the u + d 
entropy. Taking this at a face value changes the result 
(||) by less than 10 %. 

In conclusion, we showed in this paper that detection 
of QGP formation is quite possible through the simple 
measurement of h + /h~ fluctuation. This measurement 
should be very feasible for STAR. We also emphasize that 
this is a direct confirmation of the lattice QCD results. 

What we considered here is the simplest ratio out of 
many possible ones that can behave quite differently in 
the presence of a QGP. For instance, the strangeness 
anti-strangeness ratio fluctuations can provide us with 
a valuable handle on the strangeness distributions with 
or without the formation of a QGP. These and related 
issues are under active investigation and will be reported 
elsewhere 0. 
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Note added: After finishing this work, we received a 
preprint by Asakawa, Heinz and Muller [ fl2| which ad- 
dresses similar issues. 
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